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ABSTRACT 

We present the spatially-resolved observations of HCN J = 1 - 0 emission in 
the nearby spiral galaxy M51 using the IRAM 30 m telescope. The HCN map 
covers an extent of 4' x 5' with spatial resolution of 28", which is, so far, the 
largest in M51. There is a correlation between infrared emission (star formation 
rate indicator) and HCN (1-0) emission (dense gas tracer) at kpc scale in M51, 
a natural extension of the proportionality between the star formation rate (SFR) 
and the dense gas mass established globally in galaxies. Within M51, the relation 
appears to be sub-linear (with a slope of 0.74±0.16) as rises less quickly 
than Lhcn- We attribute this to a difference between center and outer disk 
such that the central regions have stronger HCN (1-0) emission per unit star 
formation. The IR-HCN correlation in M51 is further compared with global one 
from Milky Way to higli-z galaxies and bridges the gap between giant molecular 
clouds (GMCs) and galaxies. Like the centers of nearby galaxies, the L IR /L H cn 
ratio measured in M51 (particularly in the central regions), is slightly lower than 
what is measured globally in galaxies, yet is still within the scatter. This implies 
that though the Lir/Lhcn ratio varies as a function of physical environment in 
the different positions of M51, IR and HCN indeed show a linear correlation over 
10 orders of magnitude. 

Subject headings: galaxies: ISM — radio lines: galaxies — ISM: molecules — galaxies: 
individual(M51, NGC 5194) 
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INTRODUCTION 


The molecular gas in galaxies is organized in giant molecular clouds (GMCs). Within 
these clouds, denser regions, often called chimps and cores, are found and that is where 


the star formation takes place. Th e star 


the amount of molecular gas flKennicutt 


orma 


1989 


ion rate (S 


Bigiel et ah 


R) is linearly correlated with 


20081 ) with the exception of 


the infrared luminous galaxies, whose SFR is higher than would be expected from the CO 
emission (molecula r gas tra cer). The SFR and total gas (H 2 and HI) surface densities are 


correlated flKennicutt 


19981) . typically with power-law exponents greater than one when 


the HI dominates, i.e. the lower column de nsity regions generally in the outer disk, and 


approaching 1 in the H 2 dominated regions (IBuat et ah 


1989 


Buat 


1992 


Bigiel et ah 


200811 . 


Thus, at the whole-galaxy scale, the important link is between the SFR and molecular gas 


content, particularly the mass of dense gas. The star 


is hig 


rer in galaxies with a higher dense gas fraction flGao fe Solomon 


'ormatioii efficiency (S 


2004a). 


E; SFR/CO 


Lada et ah 


( 120 121) suggest that the SFR and molecular gas correlation is always linear with the same 


dense gas fraction from clouds to galaxies. 

A few environments exist in which the dense gas fraction, typically the HCN/CO 
i ntensity rati o, is_clearly different from the average value of normal galaxies (~ 1/30, 


Gao fe Solomon 


2004a ). Ultraluminous infrared galaxies have a higher dense gas fraction, 


while environments having un dergo n e a gas-gas co llisio n typi cally have a lower dense gas 
fraction (e.g. 


Gao & Solomon 

2004a: 

Braine et al. 

2003, 

2004) 


The physical conditions of active star-forming GMC cores are revealed by emission 
from higli-J (>3) rotational tr ansit ion CO lines and large dipole-moment molecules like 


CS, HCO + , and HCN (Evans 1999 k The HCN J = 1 — 0 transition traces denser gas 


(> 3 x 10 4 cm ~ 3 ) than CO J = 1 - 0 (> 300 
From a HCN emission survey of 65 galaxies, 


cm 3 ) because o 


its_larger dipole moment. 


Gao & Solomon ( 2004al ) find a tight linear 
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correlation between the SFR, traced by infrared (IR) luminosity (Fir), and dense gas mass 
(Mjense), traced by the HCN luminosity (Lrcn), with an almost constant average ratio 
Fjr/Fhcn = 900 L^/( K km s -1 pc 1 2 ), whereas the IR-CO correlation is clearly non-linear. 


Wu et ah 


(12005 


20101) extended this correlation to Galactic GMCs in the^hlky Way. The 


HCN-IR correlation appears to be valid at high redshift as well ( Gao et al 


l arge gap is present between global galaxies and Galactic GMCs (cf. Figure 2 in 
20121 ) . 


2007), but a 


Lada et al. 


Few HCN m aps of externa l ga 


Gao fe Solomon 


2004b 


Kepley et al 


axies are available (e.g., 


Nguyen et al 


1992 


2 0141 ) because the HCN emission in normal spiral 


disks is weak, such that little is known about the dense gas distribution outside th e cent ral 


regions in external galaxies. Except for a few hig 


l-resolution HCb 


1996; 

Heifer & Blitz 

1997; 

Scliin 

lerer et al. 

2007; 

Krios et al. 

2007; 

Muraoka et al. 

2009; 

Lerov et al. 

2015 

) of the central kpc and a few 


maps ( 

Kohno et al. 

Levine et al. 

2008; 


342 (40" x 60") and M83 (40" x 50") and M51 (1'.5 x 2', 


Nguyen et al¬ 


e-dis h maps of IC 


19921) . few HCN 


observations exist toward the outer disks. More observations of HCN emission in nearby 
galactic disks are necessary in order to properly interpret whole-galaxy data which is a 
mixture of nucleus and disk emission. 


M51 is an ideal target to study the spatially resolved IR- 


4CN correlation across the 


It is a large local spiral galaxy ID ~ 7.6 Mpc 

(i ~ 22°, 

Colombo et al. 

2014 

). The IlerschejJ 


Ciardullo et al. 


2002) seen 


Nearby Galaxy Su r vey 


(VNG S) mapped M51 at wavelengths of 70, 160, 250 and 350 um ( Mentuch Coo per et al. 


2 012). CO emission 


Schuster et al. 


ras been mapp ed over the entire galaxy by 


(12007 


Koda et al 


CO 2-1) and 


Vlahakis et al. 


(2011 


, CO 1-0), 


(120131 . CO 3-2). The PdBI Arcsecond 


1 Herschel is an ESA space observatory with science instruments provided by European-led 


Principal Investigator consortia and with important participation from NASA. 
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Whirlpool Survey (PAWS, 


Schinnerer et al. 


20131 ) observed the CO 1-0 emission in the 


central 9 kpc of M51 at 1" resol ution. HCN emission in the central disk (1'.5 x 2') has been 


mapped by 


Nguyen et al 


M51 are explored by 


(1992b The correlations between S 


Kennicutt et al. 


(1200711 and 


Bigiel et al. 


R and H 2 at kpc scale in 


( 120081 1 . 


In this paper, we present a map of M51 in the HCN J = 1 — 0 line at 88.6 GHz with 
the IRAM 30 meter telescope to obtain a spatially resolved HCN image of a nearby spiral 
galaxy covering 4' x 5'. The IRAM 30m beam corresponds to a spatial resolution of ~1 
kiloparsec and samples large GMC associations and individual GMCs within the beam 
along the spiral arms to starburst-like regions near the center. Combined with the Herschel 
VNGS 70 /mi and 160 /im maps, we can make a detailed point-by-point comparison of the 
IR-HCN correlation at kpc scales within a normal star-forming disk, bridging the luminosity 
gap between galaxies and individual GMCs. 


The Observations and data reduction are presented in Section 2. The HCN maps and 
the kpc-scale IR-HCN correlation are presented in Section 3. In Section 4, we compare the 
distribution of HCN, CO and IR, and discuss the effect of dust temperature on the HCN 
emission. The IR-HCN correlation in M51 is compared with existing data for galaxies and 
Galactic clouds. 


2. OBSERVATIONS AND DATA REDUCTION 

2.1. HCN Observations 

During May 2004, June 2005 and July 2007, we used the IRAM 30 meter telescope on 
Pico Veleta near Granada, Spain, to map the HCN J = 1 — 0 emission, at a rest frequency 
v = 88.63185 GHz, from M51. The angular resolution of the IRAM 30m is 28" at 88.6 GHz, 
which is ~1 kpc at the distance of M51. The map covers 4' x 5' (8 kpc x 10 kpc) and the 
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spectra are spaced by 15" (see Figure 1). So far, this is the largest HCN map of M51. 


Pointing were checked about every 2 hours on nearby QSOs with strong millimeter 
continuum emission or occasionally planets, which provided a pointing accuracy of better 
than 3". The A100 and B100 receivers which receive two orthogonal polarizations were used 
to observe the molecular lines simultaneously. We used the 512 x 1MHz filter bank backend, 
yielding a spectral resolution of 3.3 km/s. Typical system temperatures at 88.6 GHz were 
~130 K on average on the T\ scale. To convert the observed antenna temperatures to the 
main-beam temperature scale, we use T m b = Xa* F eS /B cS where the forward efficiency F eS 
= 95% and the beam efficiency B e g = 78%. Observations focused ou the spiral arms for 
the outer disk regions due to the limited observing time available. Some bad weather data 
were discarded and the total observing time of the useful data in M51 is ~ 81 hours. The 
integration time at each position was 6-90 minutes depending on the signal strength. 


All the data were reduced with the CLASS program of the GILDAS package!. After 
eliminating bad channels, a linear baseline is fit to correct the total power variations in 
the receiver and atmosphere. A small fraction of the spectra are not well fitted by a linear 
baseline. The line windows are determined from the HCN spectra themselves when the 
emission is strong. Where the emission is weak or not detected, we use line windows based 
on the CO emission, which is detected throughout the map. Individual spectra with poor 
baselines are discarded and spectra at each position are combined. 


2.2. CO Data 


The CO J = 1 - 0 total power and short spacing data were taken with the 5x5 Beam 


Array Receiver System on the Nobeyama Radio Observatory 45 m telescope fiKoda et al 


2 http://www.iram.fr/IRAMFR/GILDAS/ 
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2Dllh . 

All CO and the Herschel IR data (70, 160 and 250 //in) are convolved to the resolution 
of the HCN observations (28") to make a detailed point-by-point comparison with HCN. 


3. RESULTS AND ANALYSIS 

3.1. HCN line 

3.1.1. HCN Spectra 


The HCN J = 1 — 0 spectra are shown in Figure 1. All spectra are plotted in main 
beam temperature (T m b) and smoothed to 13.3 km resolution for show (for some weak 
HCN positions, the spectra are further binned to 26.7 km s” 1 for display). Figure 2 shows 
several HCN and CO line profiles in detail. The HCN emission line is centered at 475 
km s” 1 in the nucleus and shifted to 535 and 400 km s -1 in the northern and southern 
regions respectively. The peak of HCN spectrum is 1/10 of the CO peak in the center 
regions while it is 1/30 in the outer disk. The emission windows of the HCN lines are 
very similar to that of CO lines and the line widths tend to decrease as the galactocentric 
distance increase. In the central regions, the zero power line width is close to 200 km s -1 
and decreases to 30 — 60 km s” 1 in the outer disk. The HCN li ne peak an d profile in the 


center are consistent with the previous observations of 


Nguven et ah 


(119921 ) 


3.1.2. Local HCN Line Luminosities 

The HCN integrated intensities were measured as /hcn = f T mh dV over the emission 
window. For the HCN spectra with high signal to noise ratio (S/N > 5) in integrated 
intensity, we define the emission window from HCN. But for the HCN spectra with low 







100 0 -100 

Aa(arcsec) 


Fig. 1.— The HCN J = 1 — 0 line profiles observed in M51. Each profile is represented 
as a function of the antenna pointing offset in right ascension (RA) and declination (DEC) 
relative to the map center (a = 13 h 29 m 52.7 s , S = 47°11 , 43.0"). The vertical and horizontal 
axes of each spectrum are T m b (K) from -0.01 to 0.06 mK and V (km s _1 ) from 320 to 620 
km s -1 , respectively. All the data have been smoothed to a velocity resolution of 13.3 km s^ 1 
for show (for some weak HCN positions, the spectra are further binned to 26.7 km s -1 for 
display). 
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Fig. 2.— Spectra of HCN (black solid lines) and CO (blue dashed lines) at representative 
positions in M51. All spectra are on the T m b scale and binned to a velocity resolution of 6.7 
km s' 1 for HCN and 5.2 km s” 1 for CO. The temperature scale of CO spectra are divided by 
30 for comparison purposes. Coordinates of the observed positions (An, A5) are indicated 
at the top right and the HCN integrated line windows are shown at the bottom in each box. 
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intensity S/N (< 5), we define the emission window using the CO spectra, detected in all 
positions. As the intensity is depend on the strength of line relative to zero point d efined 


by the baseline outside the line window flBraine et a l 


1993 


Matthews fe Gao 


20011) . the 


integrated intensity errors were estimated as 


^hcn — r rms \/Whcn8c/(1- ~ Whcn /W), 


where T rms is the rms noise fluctuation, Whcn is the width of the line window in km s -1 , 
S c is the channel spacing in units of km s _1 , and W is the entire velocity coverage (1000 
km s _1 ). 3<5hcn upper limits were established for the weak S/N < 3 spectra. 

The HCN luminosity in each observed region is calculated as: 


^hcn(A" km s 1 pc 2 ) = 23.5 kld 2 L (Mpc) /hcn(A' km s 1 ), 


wh ere 11 is the solid angle of the regions covered by the telescope beams in arcsec 2 


([Solomon et al. 


19971 ) and d l is the luminosity distance. The local CO luminosity is 


calculated with the same formula (switching HCN with CO). 


3.1.3. HCN Distribution 

Figure 3 shows the CO, HCN and 70 /irn integrated intensity maps respectively. 

The HCN integrated intensity peak is 5.54±0.04 K km s _1 , located at 15" west of the 
optical nucleus and about 15" north of the CO peak (relative to the HCN resolution). The 
HCN map reproduces the large-scale structure traced by CO fairly well. The northern 
and southern CO arms are detected in HCN. 70 /irn image (traces SFR) shows the same 
structure with HCN and CO. The high resolution 70 /irn peaks in spiral arms are tightly 
relate to HCN peaks. 
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Fig. 3.— (a) Map of CO 1-0 integrated intensity [K km s x ] showing M51 fjKoda et al. 
2 011 ): contours: J co — 14-7 24.7 34.7 44.7 54.7 K km s _1 ; top: color scale in K km s -1 


(b) Map of HCN 1-0 integrated intensity [K km s -1 ]; contours: /hcn = 0.1, 0.6, 1.9, 3.4, 
4.9, 5.4 K km s -1 ; top: color scale in K km s -1 . (C) Map of 70 /.irn intensity [Jy/pixel] 
(with a pixel size of 1.4"); contours: / 70 /mi = 0.003, 0.009, 0.027, 0.081 Jy/pixel; top: color 
scale in Jy/pixel. The beam size is shown in the lower left corner. The cross marks the [0,0] 
center position at a = 13 /l 29 m 52.7 s , 5 = 47°11 , 43.0" (J2000). All the antennae temperatures 
involved in the integrated intensities are on the T m b scale. 

















iog(L IR )[y 
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Fig. 4.— IR (y-axis) as a function of HCN 1-0 and CO 1-0 (x-axis). Each point corresponds 
to an individual pointing. The central regions (filled triangles) and outer disk (open triangles) 
regions are distinguished in these correlations. The solid lines show the best fit. The dashed 
line sh ows consist ent Ltr/Thon value (900 K km s^ 1 pc 2 ) derived from the global galaxies 


flGao fe Solomon! 


2004a). Dotted lines indicate average ratio of IR/HCN in the outer disk. 
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3.2. Correlation between HCN and IR 

3.2.1. Lir from 70 /im and 160 pm fluxes 


In order to compare the IR and HCN point by point, we need to calculate the IR 
luminosity of each HCN beams. The Herschel PACS observations at 70 and 160 //in bracket 
the peak of the thermal dust em ission and a s such provide a n excel lent estimate of the dust 


luminosity (Boq uien et ah 


20111). Following 


Galametz et al. 


(120131) . we calculate the total 


IR luminosity from the convolved 70 //m and 160 //m images as: 


log Lir = [(0.973 x log zaLj,( 70) + 0.567) + (1.024 x log z/L v (160) + 0.176)]/2. 


3.2.2. Correlation between HCN and IR 

Figure 4 shows Lir vs. Thcn f° r all of the points with HCN detection ( S/N > 3). 
The HCN emission is correlated tightly with the IR over 1.5 orders of magnitude in HCN 
intensities (correlation coefficient R = 0.87). The least-square fits yield a slope of 0.74, such 
that the SFR varies less than the dense gas mass. The best-fit (logarithmic) relation is: 

log(L IR (L Q )) = (0.74 ± 0.16) log(L H CN(A' km s -1 pc 2 )) + (4.28 ± 0.95). (1) 

The uncertainties given in the equation refer to the fitting errors responding to the 
uncertainty of Thcn- The dense gas star formation efficiency (Lir/Lhcn) of the central 
regions (filled triangles, with a average Tir/Thcn ratio of 388 ± 47 L e (K km s" 1 pc 2 ) _1 , 
just the HCN luminosity errors are taken in account) is systematically lower than that 
of the outer disk regions (open triangles, with a averaged Tir/Thcn ratio of 691 ± 156 
L q (K km s -1 pc 2 )" 1 ). 
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3.2.3. Comparison of the HCN-IR and CO-IR correlations 

The relation between CO and IR luminosity of the same regions in M51 is also shown in 
Figure 4. The least-square fit yields a slope of 1.07±0.28 (correlation coefficient R = 0.78), 
such that the HCN varies more than the IR and CO for the same regions. The correlation 
coefficient between HCN and IR (0.87) is better than that of the CO-IR correlation (0.78). 
The fitting error for the slope of IR-CO (0.28) is much larger than that of IR-HCN (0.16). 
The most obvious difference occurs at the high Lir end, dominated by the central region, 
similar to the situation for galaxies as a whole. In Figure 4, lines of constant IR/HCN and 
IR/CO ratios fit the outer disk regions, but the constant IR/HCN line lies above most of 
the central regions, while the constant IR/CO line does not. 


4. DISCUSSION 

4.1. Comparison of the HCN, CO and IR Distribution 

The azimuthal averages are calculated to investigate the radius distribution of CO, 
HCN and IR displayed as lines in Figure 5. We average over rings of 15" width (regions 
within the intervals r = n x 15" + 7".5 ± 7".5, n = 0,1...10), comparable to the angular 
resolution of our data. All the intensities of CO, HCN and IR decrease rapidly with radius 
in the central 45"(1.66 kpc) and then more slowly at larger radii. HCN decreases faster 
than IR, and IR decreases faster than CO along with radius. The HCN emission is more 
confined to the inner regions than IR and CO. 

The HCN/CO, IR/HCN (SFE^ pns p ) and IR/CO (SFE) distributions are shown in 
Figure 6. The HCN/CO and IR/CO ratios are much higher in the center than in the outer 
disk, but IR/HCN is lower. The IR/HCN ratio is 388 L@(K km s” 1 pc 2 ) -1 on average in 
the central 45"(1.66 kpc), while it is 691 L@(K km s -1 pc 2 )” 1 in the outer disk, illustrating 
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R [arcsecl 

0 50 100 150 



Fig. 5.— Intensity of HCN in the central regions (filled triangles) and in the outer disk (open 
triangles) as a function of the offset from the assumed center of M51 (a = 13 ft 29 m 52.7 s , 6 = 
47°11 , 43.0"). Average values of CO (blue, divided by 30), HCN (black) and IR (red) within 
the intervals r = nx 15" + 7.5"±7.5" where n = 0, 1... 10 are displayed by lines. Dotted line 
show the radii of 45"(1.66 kpc) to distinguish the center and outer disk. 
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Fig. 6.— (a) HCN J = 1 - 0 integrated intensity map [K km s -1 ] of the detected regions 
(S/N > 3). (b) (c) (d) are respectively maps of L H cn/-^co, -^ir/^hcn [Lq/(K km s -1 pc 2 )} 
and Li R /Lco [L q /(K km s" 1 pc 2 )} of the same region. Contours are the same as those in 
Figure 3(a). The black circle shows the assumed center occupying the central 45"(1.66 kpc). 
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that the physical conditions in the center are quite different from the outer disk. 

A pattern of high L m /L C o (SFE) and L ir /L H cn (SFE dense ) can be seen in some regions 


of the outer spiral arms (Figure 6 

. GMAs cores in the spiral arm are found to have nearbv 

star-forming regions ( 

Eenif 

a et al. 

2011 

) which suggests that star formation can be triggered 

in spiral arm. 

Koda et al. 

(2012 

) study CO lines in M51 and find that the CO(2-1)/CO(1-0) 


ratio is enhanced in the spiral arms, especially downstream of the molecular arms. Higher 
resolution and sensitivity HCN observations are needed to study the effect of spiral arms on 
the HCN-IR correlation. 


4.2. Mass of Dense Molecular Gas Traced by HCN 


As the measured IR-HCN correlation physically comes from the connection between 
star formation activity and dense molecular gas. It should be noted that the mass of 
molecular gas estimated from HCN include considerable uncertainties because HCN 
emission is optically thick. Sandstrom et al. (2013) find that in the HERACLES sample the 
CO emission to H 2 factor (M H2 /L co ) increases from t. 
out. in the disk, in agree ment with earlier work (e.g. 


re central kiloparsec 

;o regions further 

Nakai & Kuno 

1995; 

Sodroski et al. 



1995 


Brainc et al. 


19971 ). This suggests that M dense /LucTss may increase with galactocentric 


radius as well. In this situation, the SFR-M dense correlation may be linear while L TR — L HCN 
may not, requiring a varying M dense /L H cn factor to convert from L H cn to M dense . 


Gao fe Solomon! fl2004bl ) suggest that M denS e/L H cN may vary linearly with the gas 


excitation temperature. They estimate that for an excitation temperature of 35 K, 
M d ense/AHCN ~ 10 M & (K km s- 1 pc 2 ) -1 . The position-dependent Mdense/ Lhcn conversion 
can then be estimated as: M dense /-LHCN = 10 x (35 /T)M e (K km s _1 pc 2 )~L Since we have 
no means of estimating the excitation temperature, we assume that the dust temperature is 
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representative of the temperature of the gas in dense regions. This is probably not entirely 
unreasonable as the dust and gas exchange energy fairly quickly at high density. The dust 
temperatures are estimated from the Herschel 160 /mi and 250 //m fluxes assuming an 
emissivity index of /3 = 2. While most spirals show a clear decrease in dust temperature 
with radius, in M51 we do not find a radii trend and the dust temperature varies little 
(24-28 K) in the central disk and outer spiral arm regions, resulting in little change in 

Mdense / -^HCN- 


There are still some weaknesses for our estimation of M dense /L H CN- The gas excitation 
temperature may be different from the dust temperature in GMC cores. And we do not 
consider the decreasing of HCN abundance which may be the primary cause of the apparent 
increase of M dense /L H cn with radius. Another possible cause for the higher average IR/HCN 
ratio in the outer disk is that some of the IR emission may come from the general interstellar 


radiation field, w 


rich probab ly cont ributes proportionally more to the outer disk IR than 


to the center IR (ICalzett i et al . 


200511 . 


4.3. Spatially-Resolved Star Formation Laws 


The spatially-resolved star formation laws quantify the relationship between the 
SFR and molecular gas in nearby galaxies at scales of ~ 1 kpc or less . SF R-H 2 relation 


shows various power-law in dexes, ranging 


et al. 2008), ~ 1.1 - 1.5 ( Kennicutt et al 


fom < 1 (e.g. 


2007; 


Shettv et ah 


Schruba et al 


20131 1. - 1 (Bigiel 


20101 1. A single linear 


spatially-resolved SFR-H 2 relationship can not describe relation between star formation 


and molecular gas for the large scatter from galaxy to ga laxy fjSaintonge et ah 


2011 


Schruba et al. 

2011; 

Rahman et al. 

2012; 

Lerov et al. 

2013) 


Correlations between SFR and dense molecular gas tracers, such as HCN 1-0 
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(Liu. Gao fe Greve 


(Liu et al 


2015|), HCN 4-3, HCO+ 4-3, CS 7-6 jzhang et al.l l2014h and high-J CO 


20151) . show a tight linear correlation in global galaxies. However, IR-HCN shows 


a sub-linear relation in M51 at kpc scale which is likely biased by the limited dynamical 
range, particularly the high IR central regions. And it is possible that the galaxy variations 
may dominate this inconsistency as shown in SFR-H 2 relation. More HCN survey toward 
galaxy is needed to explore the spatially-resolved IR-HCN relation. 


4.4. L 


HCN 


vs. L ir : Correlation from GMCs to High-z Galaxies 


One of the motivations for studying HCN in M51 comes from the intimate link 
between star formation and dense gas in galaxies and in individual GMC cores in the 
Milky Way. I 11 Figure 7, we show plots of Lir versus Lrcn and Lir/Lrcn (an indicator 
of SFE dense ) versus L IR measured locally in M51 (center, filled black triangles; outer disk, 


open triangles) 


compared with the observations of individ ual GMC cores in t 


(open squares; 
(£lled__blue squar es; 


Wu et ah 


2010),M33 (hi 


Brouillet et ah 


ed black squares; 


re M ilky Way 


Buchbender et al.l 


20131 ). M31 


Krips et al 


20051) . nearby g;alactic centers (filled b lue triangles; 


20081), whole 


galaxies (filled hexagons; 


galaxi es ( open 


Gao ct al. 


hexagons; 


Gao & Solomon 


2004bf ) and high-z 


2007 ). For whole galaxies, the mean value of Lir/L H cn 


is 9OOL 0 (A” km s -1 pc 2 ) 1 (shown as the solid line), while it is 629L 0 (A' km s” 1 pc 2 )^ 1 
in M51. The Lir/Lrcn ratio in M51 is a factor 1.5 lower on average than that of most 
galaxies, yet is still contained in the scatter. This suggests that a linear correlation is 
highly applicable to different physical scales, linking GMCs to kpc local regions and whole 
galaxies in a single relation. The linearity of the HCN-IR correlation was interpreted as a 
fundamental unit of all star-formin g syste m, thus both IR and HCN are simply accumulated 


by adding in more units flWu et al 


20051) 


I 11 M33, 


Buchbender et al 


(120131) observed 7 clouds at roughly lOOpc resolution in 
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CO, HCN and HCO + . It is apparent from Figure 7 that the HCN emission for a given IR 
luminosity is weak in M33 (filled black squares), a factor 3 on average weaker than in the 
disk of M51. A reason to expect the HCN line to be weak with respect to IR is the sub-solar 
metallicity in M33. The dust (where IR from) abundance is expected to vary approximately 
linearly with metallicity (i.e. O/H) assuming a similar fraction of heavy elements in dust 
phase. Because N abundance decreases more quickly than the O abunda nce ( i.e. N/O 


ratio decreases with decreasing O/H) (see Figure 18 of review by 


Garnett 


20021 ). the HCN 


line is expected to decrease more than IR relative to Galactic abun dances in M33. The 
disk of M31 was obs erved at a s imil ar r esolution by 
luminosity from the 


Brouillet et ah 


Draine ct ah 


020051 ): using the IR 


020141) maps, these points are also shown in Figure 7 


(filled blue squares) and are in rough agreement with the M51 IR/HCN ratio. From the 
resolved observations, the IR/HCN ratio is low for galactic centers, moderate (i.e. following 


the Gao et ah relation) for the disks of M51 and M31, and high fo r M3 3. T 


is covered by the individual galactic clouds studied by 


Wu et al 


(2010 


2005 


lis whole range 
)• 


The centers of 9 nearby galaxies NGC 1068, 2146, 3627, 4569, 4826, 5194, 6946, M82 


and Arp 220 have been observed in HCN(J = 1-0) with the IRAM 30 m ( Krips et al.l 20080 . 


These objects exhibit nuclear star burst and/or AGN signatures. We estimate the infrared 


luminosity of the centers from the convolved Herschel 100 /mi map flGalametz et ah 




20131 ) 


observed by KINGFISH_| or the combination of the convolved Herschel 70 and 160 //m map 
(see section 3.2.1) observed by VNGS. NGC 2146 and M82 centers show higher Li R /L H cn 
ratios than M51, the other 7 nearby galactic centers are consistent with the central regions 
of M51. The lower Ti R /Thcn ratios of nearby galactic centers compare to the mean value 
of whole galaxies (combination of outer disk and center ) suggest that the L ir /L H cn ratio 
may be higher in the outer disk than the center, which support the increasing trend from 


3 Key Insights on Nearby Galaxies: a Far-Infrared Survey with Herschel 
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center to outer disk in M51. More HCN observations toward the outer disk are needed to 
study this trend. 


4.5. Li R /Lco vs. Rhcn/T 


co 


Figure 8 shows a correlation between Lir/Lco (tracer of star formation efficiency; SFE) 
and Lhcn/Fco (tracer of dense gas fraction; -Rhcn/co)- The SFE in the central regions is 


syste matically higher than the outer disk regions, which agrees well with 


Gao fe Solomon 


(I2004ai ) that the star formation efficiency depends on the fraction of molecular gas in a 


dense phase. However, the correlation between SFE and -Rhcn/co is not linear. -Rhcn/co 
in the center is higher than the linear fit determined from the mean value of the outer disk 
regions. Thus, either (i) Lir underestimates the SFR in the center (relative to the outer 
disk) or (ii) L H cn overestimates M dense in the center (relative to the outer disk) or (Hi) the 
central environment makes star formation less efficient. The 24 pm emission is less sensitive 
to the general interstellar radiate held and shows the same trends as in Figure 4, so we do 
not support (f). Our analysis excludes temperature as the primary cause of (ii). Abundance 
variations, presumably driven by Nitrogen, could provide at least a partial explanation for 
the trends we fold. Our work cannot address (in). To further study the effect of density 
on star formation, we are actually observing the higher dense gas tracer (HCN 4-3) in the 
central disk of M51 with the JCMT telescope, p otentially extended the IR and HCN 4-3 


relation in the global galaxies ( jZhang et al. 


20141). 


SUMMARY 


We present the spatially-resolved observations of the HCN J = 1 - 0 transition in the 
nearby spiral galaxy M51 using the IRAM 30 m telescope. The map covers an extent of 
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Fig. 7.— Fir-Lhcn correlation (upper) and Fir/Fhcn-Fir correlation (lower) for M51 
observations (open and filled triangles), M33 (filled black squares), M31 (filled blue squares), 
nearby galactic centers (filled blue triangles), whole galaxies (open hexagons), liigh-z galaxies 
(filled hexagons) and cloud cores (open squares). 
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I°9(L|-icn/Lco) 


Fig. 8.— Correlation between Lir/Lco (trace SFE) and Lrcn/Fco (-Rhcn/co) for the re¬ 
gions of M51 and global galaxies. Central (filled triangles) and outer disk (open triangles) re¬ 
gions in M51, galaxies (open blue hexagons) are distinguished in this figure. Solid li ne shows 


th e cor relation between Lir/Lco and LncN/-bco for the global galaxies from 


Gao fe Solomon 


( 2004a ) and dotted line show the consistent ratio of Lir/Lco and LncN/-hco for the outer 


disk regions. 
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4' x 5' with spatial resolution of 28" and 15" sampling, which is, so far, the largest HCN 
map of M51. The HCN emission peaks 15" west of the optical nucleus and about 15" north 
of the CO peak (at the HCN resolution). 

Comparing the individual observed positions, the IR emission varies less than the HCN 
emission, such that Lm oc L^^ 0 ' 16 , despite the overall linear relation between HCN and 
IR emission from Galactic clumps to distant galaxies. While the conversion from HCN 
intensity to dense gas mass may vary, temperature (as traced by the dust) cannot explain 
the variation. It is possible that the galaxy variations may dominate this inconsistency as 
shown in SFR-H 2 relation. 

Lin traces SFR and Lhcn traces dense gas mass, so the ratio provides information 
about the efficiency of transformation from dense gas into stars. There is a clear radial 
trend in L ir /L H cn decreasing with galactocentric radius. Comparing with resolved data in 
other galaxies, 7 nearby galactic centers show a consistent Tir/Thcn ratio (an indicator of 
SFEdense) with the central regions in M51. While the ratio measured locally in M51 is 1.5 
times lower than the galactic average, it is still within the scatter. Though Tir/Thcn may 
vary with physical environments in M51, HCN and IR indeed shows a linear correlation 
over 10 orders of magnitude. 


We are very grateful to Jin Koda for providing the CO 1-0 data on NGC 5194. The 
VNGS data was accessed through the Herschcl Database in Marseille (HeDaMcl) operated 
by CeSAM and hosted by the Laboratoire d’Astrophysique de Marseille. We appreciate the 
generous support from IRAM staff during the observations and data reduction. This work 
is partially supported by NSFC grant Nos. 11173059, 11390373, 11273015 and 11133001), 
CAS No. XDB09000000 and the National Basic Research Program (973 program No. 
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2013CB834905). 
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